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Reaction dynamics involving subdiffusive species is an interesting topic with only few known re- 
sults, especially when the motion of different species is characterized by different anomalous diffusion 
exponents. Here we study the reaction dynamics of a (sub) diffusive particle surrounded by a sea of 
(sub)diffusive traps in one dimension. Under some reasonable assumptions we find rigorous results 
for the asymptotic survival probability of the particle in most cases, but have not succeeded in doing 
so for a particle that diffuses normally while the anomalous diffusion exponent of the traps is smaller 
than 2/3. 



I. INTRODUCTION 



^•^ _ In the traditional version of the trapping problem, a normal diffusive (Brownian) particle (A) wanders in a medium 

r^ . doped at random with static traps (B), and disappears when they meet, A + B ^ B. The quantity of interest is 

c/j ' the survival probability P{t) of the particle A. This problem dates back to Smoluchowski's theory of reaction rates 

■4—^ . at the beginning of last century, and is one of the most widely investigated and applied problems of non-equilibrium 

3 ' statistical mechanics jH Q, y, |j, IM HI - An important variation of the basic trapping problem, in which a diffusive 

^ . particle wanders in a medium in which the traps are also diffusive, has been the subject of intense research since the 
seminal work of Toussaint and Wilczek [3 ■ 

The principal quantity of interest in the trapping problem is the survival probability P(t) of the A particles. 

Q From this survival probability one is able to calculate essentially all other quantities of practical interest. Yet this 

CJ ' probability is usually difficult to calculate, and the few instances in which it has been obtained are considered landmark 

'~~', contributions. In 1988, Bramson and Lebowitz |3,|3| proved rigorously that the long-time survival probability P{t) of 

I ' a particle diffusing in a one-dimensional medium doped with diffusive traps decays as P{t) ~ exp(— Ai^"), A being an 

^ , undetermined parameter. The evaluation of this constant proved elusive for many years, engendering much confusion 

00 ■ and proposed solutions that were mutually contradictory. Finally, quite recently Bray and Blythe |l(J . \l^ proved in 

^^ [ a simple and elegant way, assuming the so-called "Pascal principle," that the survival probability P{t) of a diffusing 

\l ■ particle with diffusion coefficient D' in a d-dimensional medium with d < 2 in which the traps are also diffusive 

^T I with diffusion coefficient D is independent of D' for long times, and coincides with the survival probability of an 

St . immobile target {D' = 0) in the presence of a density p of diffusive traps. In particular, in a one-dimensional medium 

r-^ ■ P{t) ~ exp(— 4p(Dt/7r)-'^/^). Bray and Blythe obtained their results by calculating an upper and a lower bound for 

— «. , the survival probability that converge to one another asymptotically. Some (but not all) of the bounding results of 

j^ ■ Bray and Blythe have been extended by Oshanin et al. [12i| to systems where the traps perform a compact exploration 

Ch I of the space, i. e., where the fractal dimension d^ of the trajectories of the traps is greater than the dimension d of 

I ■ the space. 

' O ' The purpose of this work is to extend the procedure and results of Bray and Blythe [l^, [ij, which are valid 

5 ' for a Brownian diffusive particle and Brownian diffusive traps, to situations in which the particle and traps exhibit 

anomalous diffusion, in particular, subdiffusion. The usual characterization of anomalous diffusion of a particle is 

through its mean squared displacement x{t) for large t: 



H ; 

CZ Here Kj is the (generalized) diffusion constant and 7 is the (anomalous) diffusion exponent. Ordinary Brownian 
diffusion (7 = 1, Ki = D) follows Fick's second law, ( x^J tS) ex t. The process is called sudiffusive when < 7 < 1. 
Subiffusive process are ubiquitous in nature [isL 111, [isl lla . IitI llsL [l9l I20I |21| , and are particularly useful for 
understanding transport in complex systems J5|, l22l l23l |2J| . 

Two main approaches have been used to study subdiffusive processes. The older is based on the continuous time 
random walk (CTRW) theory with waiting-time distributions between steps that have broad long-time tails and 
consequently infinite moments, ?/'(i) ^ t~^~'' for i — > cx) with < 7 < 1. Another approach is based on the fractional 
diffusion equation, which describes the evolution of the probability density P{x, t) of finding the particle at position 
X at time t by means of the fractional partial differential equation (in one dimension) |l3l |25| 

^^P{x,t) = K,,D]-'^P{x,t) (2) 



2 

where K^ is the generahzed diffusion coefBcient that appears in Eq. Q and o -Dj ~^ is the Riemann-Liouville operator, 

Here we implement the latter approach to study the one-dimensional trapping problem in the long-time regime for 
subdiffusive (or diffusive) particles that move among a distribution of non-static traps. The traps can be either 
subdiffusive or (Brownian) diffusive. For this purpose, we generalize the ideas of Bray and Blythe pillll|. 

The problem considered in this paper is a special case of a broad class of reaction- subdiffusion processes that have 
been studied over the past decades using one or the other of the established formalisms. Using the CTRW formalism, 
Blumen et al. |2y,|23,|23 considered a variety of reactions including the trapping problem A -(- r(static) — + r(static), 
the target problem j4(static) + T -^ T, and the bimolecular reactions A + A ^ and A + S ^ 0. The moving 
particles were modeled as continuous-time random walkers with long-tailed waiting-time densities. Quite recently. 
Sung and Silbey ^29] have used the CTRW model to study the dynamics of particles that react at a boundary. A 
CTRW approach has also been applied by Seki et al. |30l y^ to study the kinetics of the recombination reaction in 
subdiffusive media. However, Seki et al. went further and, from the CTRW model, derived a fractional reaction- 
diffusion equation for the geminate recombination problem. Sung et al. |3lj| directly addressed this problem from a 
fractional diffusion equation approach, but some of their assumptions and results disagree with those of Seki et al. The 
fractional diffusion approach has recently been used to get exact solutions for two types of one-dimensional trapping 
problems: the so called one-sided problem, in which all the traps lie on one side of the particle, and the two-sided 
problem, in which the traps are located on both sides of the particle (this is the traditional or standard version of the 
trapping problem) Q^ . These recent articles |23, uM share the simplifying characteristic that the reaction takes place 
between a static particle (or fixed boundary) and a subdiffusive particle. The present paper differs from these in that 
all the reacting particles (including traps) are (sub)diffusive and, moreover, the diffusion constant and the anomalous 
diffusion exponent of each species may be different. The fractional-diffusion approach has already been employed to 
study bimolecular reactions between subdiffusive particles. In particular, the annihilation A + A ^ and coagulation 
A-\- A ^ A oi subdiffusive particles was studied |33, 13 by means of a fractional generalization of the interparticle 
distribution function method jQl . The evolution of reaction-subdiffusion fronts for A-\- B —> C reactions, where both 
A and B move subdiffusively, is also amenable to analysis by means of the fractional diffusion approach [33 . Other 
recent work on fractional diffusion and CTRW models of subdiffusive reacting particles can be found in a number of 
references HI miH H^ . 

In some cases, asymptotic anomalous diffusion behavior can be found from corresponding results for normal diffusion 
with the simple replacement of t by f wherever t appears. This can be understood from a CTRW perpective because 
the average number of jumps n made by a subdiffusive walker up to time t scales as (n) ^ f^ and, in many instances 
the number of jumps is the relevant factor that explains the behavior of the system. The simple replacement result 
is evidence of "subordination" .28,] . Examples of this phenomenon are given in Sees. 5 and 7.2 of the preceding 
reference. However, there are other instances where the behavior of subdiffusive systems cannot be found in this 
way. A simple example is the survival probability of subdiffusive particles in the trapping problem, see Sec. 5 of the 
preceding reference. In particular, for systems where each species has a different anomalous diffusion exponent, such 
a replacement becomes ambiguous. This is the case for the problem considered here. 

Bray and Blythe obtained the asymptotic survival probability of a diffusing particle in a sea of diffusing traps 
by calculating an upper and lower bound that converge asymptotically. Wc follow their procedure for subdiffusive 
particle and traps, with partial success. While it is possible to obtain convergent bounds in most anomalous diffusion 
exponent regimes, this procedure does not work in all regimes. In particular, the bounding procedure encounters 
difficulties when the particle A diffuses normally and the traps are "too slow." In Sec. ^ we calculate the upper bound 
of the survival probability, and in Sec. lIIII the lower bound. The survival probability is established, when possible, in 
Sec. II VI Section rVI and the Appendix presents a compendiary of results and some comments on open problems. 

II. UPPER BOUND FOR THE SURVIVAL PROBABILITY 

The "Pascal principle" of random walks says that the best survival strategy for a random walker A surrounded 
by a sea of trapping random walkers B is to stand still. This assumption was adopted in one dimension by Bray 
and Blythe [iflllj, and proved by Bray, Majumdar and Blythe |43| for d < 2. Almost simultaneously, Moreau et 
al. [IJ proved the Pascal principle for a rather general class of random walks on d-dimensional lattices. Although the 
case in which both particle and traps perform subdiffusive random walks was not considered, the Pascal principle is 
an intuitively plausible (and most likely provable) result for this case as well, and we will simply assume that it is 
applicable. 



The Pascal principle thus says that the survival probability Pu{t) of a static particle A surrounded by a random 
(Poisson) distribution of randomly walking traps (the "target problem" ) is an upper bound for the survival probability 
P{t) of a (sub)diffusing particle A. We proceed to calculate this upper bound for the subdiffusive target problem. This 
problem has been considered z8] by means of the CTRW model and, for the three-dimensional case |3^ by means of 
a fractional diffusion approach. Here we calculate Pu (t) by generalizing the approach of Bray and Blythe [lO, UM ^o 
the subdiffusive case. 

Consider a target of size 2L centered at the origen, and let Qi (t\y) be the probability that the trap initially placed 
at y > L has not reached the target A a,t y = L hy time t. Then J42|, in terms of the Fox's H function, 



h{t\y) = I - Hi'l 






(1,7/2) 
(0,1) 



= 1-H 



y-L 



(4) 



For 7—^1 the Fox's H function becomes the complementary error function (we set Ki — D), and the ordinary 
Brownian motion result is recovered, 



i{t\y) = 1-erfc 



y-L 

xfADi 
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1. 



(5) 



Next consider N independently diffusing traps that, at i = 0, are placed at random in the interval L < y < L + R. Here 
and henceforth 2R is the size of the system, which we will take to infinity at appropriate points in the calculations. 
The probability QN(t) that the stationary target A has survived up to time t is 



N 



rL+B. 



dy, {l-H 



Ui - L 



= <i- 



1 

R 



L+R. 



dyH 



y 



/^^ 



N 



(6) 



or, in terms of the density p = N/R of traps. 



~Xt) 



lim < 1 

R^oo 



R 



L+R 



dyH 



y/K^ 



pR 



exp <^ -p^K^ri / dzH {z\ 



We need to evaluate the integral 



10 



dz H^l 



(1,7/2) 
(0,1) 



which can be done from the properties of the Fox's PI function 43] . One finds that 

d 



H 



(1,7/2) 
(0,1) 



dz 



ttIQ 



(1 + 7/2,7/2) 
(0,1) 



But 



and 



so that 



Therefore, 



H 



(1+7/2,7/2) 
(0,1) 



10 



rriU 
-"11 



(1 + 7/2,7/2) 
(0,1) 

1 



r(i + 7/2)' 



r(i + 7/2)- 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



,{t) =exp 



r(i + 7/2) 



(13) 



This is the survival probabiHty of the target when all the traps are located to its right. When the traps are located 
on both sides of the target, the survival probability of the target is the square of Eq. H13|) : 



^c/W = QLW=exp 



r(i + 7/2) 



(14) 



This is the upper hound on the survival probability of the moving particle. 

Incidentally, as is well known, the survival probability for the target problem is related to the distinct number of 
sites S{t) visited by a trap up to time t: [23, Hj] 



Pu{t) 



^p{S{t)) 



(15) 



Comparing this expression with Eq. (|14|l . one finds that the average value {S{t)) of the territory explored up to time 
t by a subdiffusive walker with generalized diffusion coefficient Kj and anomalous diffusion exponent 7 is 



{S{t)) 



2y^K:^ 



r(l + 7/2)' 
This result agrees with that found by Yuste and Acedo [33 using a different approach. 



(16) 



III. LOWER BOUND FOR THE SURVIVAL PROBABILITY 



Let Pl (t) be the probability that the particle A (which is now allowed to move) remains inside a box of size L and 
that all the traps remain outside this box until time t. When this happens, the particle A survives. It is clear that 
Phii) is a lower bound for the survival probability P{t) of interest here because there exist many other ways in which 
the particle A may survive. The probability PL{t) is itself the product of three probabilities: 



1. The probability Qi that at t = the box of size L contains no traps: 

n, — e-pL 



(17) 



2. The probability Q2 that no traps enter the box of size L up to time t: 

2 



Q2 = exp 



r(i + 7/2) 



p^K^ti 



Puit). 



(18) 



3. The probability Q3 that the particle has not left the box of size L up to time t. We proceed to evaluate this 
quantity. 

Let W{x^t) be the probability of finding the particle A at position x at time t if it was at position a; = at time 
t = Q and there are absorbing boundaries at a; = —L/2 and x = L/2. Solving the fractional diffusion equation by 
means of separation of variables [13 one finds 



W{x,t) 



f (-1)« sin (^!i±lML±iV2)^^, ^_^,^(2^ ^ if^H^'/L- 



n=0 



(19) 



where K' and 7' are the generalized diffusion constant and the anomalous diffusion exponent of the particle A. 
Therefore, 



L/2 

W{x,t)dx 

-L/2 T^ 



^2n+l ^ 



-K',{2n+iyTrh^ /L 



(20) 



Next we distinguish two cases in the handling of the sum in Eq. (I20() : first we deal with a subdiffusive particle, and 
subsequently with an ordinary diffusive particle. In the subdiffusive case, we note that for large arguments (z 3> 1) 
the Mittag-Leffler function has the expansion 



Ey(-^)-T.m 



(-1) 



m+l 



'^ r(l - 7'm) 



(21) 



so that 



Qa 



E 



/■ -| "^rn+l r 2?n 



1 r(i — 7'™) 



1™ Z^ 



(-1)" 

(2n+l)™+i' 



Therefore, for t —f oo one finds 



(22) 



Qs 



o 



i^ 



^yi^' 



8r(i-7')if;,n' 

Consequently, a lower bound on the survival probability of the particle A is 



PUt) = Q1Q2Q3 ^ e-P"^ exp 



p^K^t-r 



L^ 



8T{l-Y)K'^,t^' 



1 + 



K'yt''' 



r(i + 7/2) 

It can easily be ascertained that this expression is maximal when L = L* = 2/p, i.e., Phit) ^ -PL*(i) with 



Pl*W = 



8r(l-7') \p) K'f 



1 



— exp 



r(i + 7/2) 



p^K-yf 



1 + 



1 



p^K'^,ty 



This is then our best lower bound for the survival probability P(t) of a subdiffusive particle. 
When the particle A diffuses normally, Eq. H2()|l becomes 



n ^-^ 2n + 1 

n=0 



with D' = K[. For long times p^[ll| 



Q3 exp [~D'nh/L^] , i > 1 



so that 



PLit)^QlQ2Q3 



-pL 



exp 



p^isr^i^ 



exp [-D'TT^t/L'^] 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 



r(I + 7/2) 
for t 3> 1. This lower bound can be again be maximized by optimizing the value of L. The optimal value is |lOllll| 



1/3 



L* = {2Tr^D't/p) ' , so that 



PLit)<PL'{t) = -exp 



2yfp^lQn 



- 3(^2^2^'t/4)'/' 



(29) 



Note that the dominant term inside the bracket depends on the value of 7, the anomalous diffusion exponent for the 
traps. We distinguish three cases: 



1. Traps with 2/3 < 7 < 1. In this case, for i ^ 1, 



2^p^K^ti 

r(i + 7/2) 



» 3(^V'^'V4)'/' 



so that 



2. Traps with 7 = 2/3. Now 



PL*(t) = -exp 

TT 



"r(l + 7/2) 



PL(0<-PL-W = -exp 



2vVi^ 

r(4/3) 



3(7r2p2i574)i/3 ti/3 



(30) 



(31) 



(32) 



that is, the second contribution in the exponent in Eq. H29|) is of the same order as the first and must thus be 
retained. 



3. Traps with < 7 < 2/3. Now the second term in the exponent of Eq. H29|) is dominant: 



Pi.(t) = -exp -3(7rV£''V4) 



\l/3 



(33) 



In the next section we examine our upper and lower bound results to establish the behavior of the survival probability 
of A whenever possible. 



IV. SURVIVAL PROBABILITY 



We now combine our upper and lower bound results. Recall that the label and exponent 7 is associated with the 
traps and 7' is associated with the particle A. The upper bound on the survival probability is in all cases given 
in Eq. H14(l . but the lower bound depends on the anomalous diffusion exponent of the particle. We distinguish the 
following cases: 

1. Subdiffusive particle (0 < 7' < 1) and diffusive or subdiffusive traps (0 < 7 < 1). The lower bound is given in 
Eq. ^, so that Fl* (0 < P(.t) < Puit) leads to 



r(i + 7/2) 



< 



< 



\nP{t) 

7^5^ -r( 1 + 7/2) 



2 In Jp^KLt-r' +2 + ln[2r(l-7') 



Vp^k:^ 



o 



P^K't^' 



(34) 



p'^K' ,tf' <C ^/p^^KI^ and the upper and lower bounds converge asymptotically. We therefore 



For t — » 00, In 

arrive at the explicit asymptotic survival probability 



P(t) - exp 



r(i + 7/2) 



p^K^ti 



(35) 



for < 7 < 1 and < 7' < 1. Note that for 7 = 1 we recover the normal diffusive result obtained earlier [lO| . 
A noteworthy result here is that the survival probability depends only on the exponent 7 that characterizes the 
traps and not on 7' that characterizes the particle. This is interesting vis a vis the subordination issue. 



2. Diffusive particle (7' = 1) and subdiffusive traps with 2/3 < 7 < 1. The bounds here are 



InF(t) 



r(i + 7/2) - ^p^ic^t^ - r(i + 7/2) 

and the asymptotic survival probability is again given by Eq. 13511 






(36) 



3. Diffusive particle (7' = 1) and subdiffusive traps with 7 = 2/3 (marginal case). Now PL'{t) < P{t) < Pu{t) 
leads to the more ambiguous inequalities 



< 



InP(t) 



< 



r(4/3) - ^^K^t^ - r(4/3) 






2/3 



j^n/3 



K, 



1/2' 



(37) 



The bounding procedure is therefore not able to predict the value of the prefactor A in P{t) = exp(— At^/^), but 
the asymptotic behavior — InP(i) c>c t^'^ is evident. 

4. Diffusive particle (7' = 1) and subdiffusive traps with < 7 < 2/3. The bounds here are also given by Eq. H36|l . 
so that the bounding procedure is not able to determine the asymptotic behavior of P{t) at all for this case. 
We are not even able to assert the asymptotic stretched exponential form P{t) ^ exp (— At'^). 



V. PANORAMA AND DISCUSSION 

Bray and Blythe [ifl UM have calculated the asymptotic survival probability of a diffusive particle A in a randomly 
distributed sea of diffusive traps B in one dimension, and have determined the precise value of the coefficient A in the 
classic result P(t) ~ exp (— Ai^/^) first obtained by Bramson and Lebowitz [alal- We have attempted to generalize 
this result to the case where one or both of the species move subdiffusively. Our particle A is characterized by the 
anomalous diffusion exponent 7' and generalized diffusion coefficient Kj>, and the traps by 7 and K^. These may be 
the first results in the literature involving two subdiffusive species with different anomalous diffusion exponents. 

When both species are subdiffusive (7 and 7' both smaller than unity), the survival probability is independent 
of 7' and determined entirely by the subdiffusive properties of the traps, cf. Eq. (|35|l . When the particle moves 
diffusively (7' = f), on the other hand, we are unable to unequivocally determine the coefficient A for all cases using 
this procedure. If the traps move sufficiently rapidly (2/3 < 7 < f ) then the result Eq. H35(l is still valid. Note that 
this reduces to the Bray and Blythe result when 7=1. The case 7 = 2/3 is marginal in the sense that we can 
establish the behavior P{t) ^ exp (— At^^'^), but are not able to determine the constant A. Note that this particular 
time dependence of the survival probability is the same as the classic result for the survival probability of a diffusive 
particle in a sea of immobile traps |43, Ha, 113, Eg . If the traps are too slow ("extremely subdiffusive"), < 7 < 2/3, it 
is no longer possible to determine even the time dependence of the survival probability on the basis of this approach. 

We thus find that in so far as one can even think of some sort of subordination principle (and whether such thinking 
is appropriate is debatable), it is determined by the behavior of the traps, i.e., by the replacement of t by f*". Even 
in the range of exponents where this is possible, it is only possible for the main asympotic contribution to P{t) but 
not for the correction terms to the leading asymptotic term. 

It is interesting to note that the value L* = 2/ p that maximizes the lower bound of the survival probability for a 
subdiffusive particle A does not grow with time. This implies that finite particle size effects could become relevant 
with increasing density p. This is completely different from the case of a Brownian particle A, since the growth 
L* ex t^''^ now suppresses such finite size contributions for any given density. 

A number of questions and opportunities for further work arise from our analysis. In our analysis (as in that of 
Bray and Blythe [iflllll) we have simply applied the Pascal principle that the best survival strategy of a particle 
in a sea of moving traps is to remain stationary. While it seems intuitively obvious that this principle would hold 
whether the motions are diffusive or subdiffusive, the proof has only been presented for the diffusive case pjl l4ll|. 
As an interesting aside, we note that while the proof of the Pascal principle has assumed an equal concentration of 
traps on either side of the particle, we conjecture that this is not a necessary condition, and that the Pascal principle 
also holds with an asymmetric distribution and even if the traps are all located on one side of the particle. This, too, 
remains to be proved. 

At this point we inject a digression that is relevant not only to our analysis but also to the original work of Bray and 
Blythe |10lllll |. They assumed that the particle A is initially surrounded by a random (Poisson) distribution of mobile 
traps, an assumption also made in our analysis, cf. Eq. (|I7|I . On the other hand, if at the start of the observations 
(t = 0) the process has already been taking place for some time — r (i.e., if the process started at some time r in 
the past), then it is known that the distribution around the surviving particles at time t = is not of Poisson form. 
Those particles that initially had close-by traps are more likely to have been trapped already than those that did not, 
so that those particles that have survived are surrounded by a region of fewer than average traps (sometimes referred 
to as a "gap"). Bramson and Lebowitz arrive at the conclusion that the configuration of B particles is nevertheless 
dominated by a Poisson random measure ^. In Appendix^ we confirm that for any finite r the gap does not affect 
the asymptotic survival probability results of Bray and Blythe. The detailed nature of the gap is different in the 
diffusive and subdiffusive cases, and unknown in the latter. However, we conjecture that it is no more pronouced in 
the subdiffusive than in the diffusive system, and that it does not affect our results either. 

Our own results of course leave a number of questions unanswered. One obvious question concerns the marginal 
role of the trap exponent 7 = 2/3 when the particle is diffusive. Why is this a marginal exponent? A connection 
between this critical value and the fact that for a Brownian particle the length that maximizes the lower bound of 
the survival probability grows as L* ~ t^^^ seems plausible, but the conceptual basis for such a relation is not clear. 

The most pressing and intriguing puzzle to resolve is that of calculating the survival probability when the particle 
A is diffusive (7' = 1) and the traps are extremely subdiffusive (0 < 7 < 2/3). Because the upper and lower bounds 
in this case do not have the same asymptotic time dependence, we are not able to say anything about this case on 
the basis of the procedures used in this paper. 



APPENDIX A: MODIFIED INITIAL CONDITION 

Suppose that the process A + B —f B with a diffusive particle and diffusive traps began at time — r, but our 
observation of the system starts at time i = 0. Even if at time — r the distribution of B's around A was random 
(Poisson), it will not be so at time t = 0. Let p{r) be the probability density of finding an empty region of length r to 
the right of A at time i = 0, so that the first trap is found to be located between r and r + dr. This function is unknown 
for our system, but one can conjecture a behavior on the basis of known (analytic and numerical) approximate results 
for the bimolecular reactions A + A ^ and A + S ^ |4flll5?il |. 



p{r) 
p{r) 



r/{r) 
exp(- 



r/{r)) 



for r/{r) < 1 
for r/{r) > 1 



(Ala) 
(Alb) 



where (r) is the mean size of the gap next to particle A at the initial time of observation. Next we need to define the 
function q{r), the probability density that a particle B is located at a distance r from A at t = 0. This function can 
be related to p{r), most easily by considering a discrete lattice and then going to the continuum limit. Let p„ be the 
probability of a gap of size n and a particle B at position n+1 to the right oi A. Thusp„ = 11" (1 "?*)%' where qi is the 
probability that there is a trap B at distance i from A. It then follows that lnp„ = ln(gi)+X]" ln(l— 9i) ~ lii(gi)— ^" qi 
and thus in the continuum limit 



p{r) — q(r) exp 



q{x)dx 



(A2) 



To find q[r) one thus needs to solve a nonlinear integral equation, a difficult task not worth the attempt given our 
imperfect knowledge of p(r). It is nevertheless straightforward to conclude that the solution q{r) = p(r) for r <C (r) 
is compatible with ljAla|l . and that q{r) =const. for r ^ (r) is compatible with (|Alb|l . For our purposes it is thus 
sufficient to approximate q(r) by 



,p7^, 0<r<(r) 



for i? — > oo. 



(A3) 



The coefficient \jR has been set by requiring that the integral of qir) from zero to R be normalized to unity. Thus, 
the presence of A affects the distribution of B's within a distance (r) from A but not beyond, where the distribution 
remains random. 

How does this deviation from a random distribution affect the survival probability? For this calculation we also 
introduce q(r; _R) = q{r)l L q{r)dr, the probability density for a B particle to be found at a distance r from A given 
that there is a i? particle somewhere in an interval of length R that begins at A. An upper bound is obtained, as 
before, by invoking the Pascal principle. Generalizing the derivation for a Poisson distribution one straightforwardly 
arrives at the expression 



Qoo{t) = lim 



lim 



q{r;R)Qi{t\r)dr 



pR 



R pR 

q{r-R)~ / q{r-R)Qi{t\r)dr 
Jo 



pR 



(A4) 



R. 



where Qi is given in Eq. jsj and Qi — 1 — Qi. Since linii^^^oo /n q{i", R)dr = 1, it follows that 



~j{t) = lim exp 

R—^OQ 



lim exp 

R^oo 



pR q{r;R)Qi{t\r)dr 



exp 



-pi? / qir)Qi{t\r)dr 



(r) ,, 



-p\ --Qi{t\r)dr+ Qiit\r)dr 



lim exp 

R — ^oo 



pR 



Jo q[r)dr Jo 



qir)Qi{t\r)dr 



(A5) 



where we have introduced the exphcit form (|A3|) in the last Hne. The integrals can be carried out in the long time 
limit, to yield 



Qoo(i) = exp 




for f -> oo. (A6) 



Since (r) is a constant (no matter how large), eventually the second term in the exponent dominates, and the 
upper bound on the survival probability in the diffusive case is exactly the same as obtained for an initially Poisson 
distribution of B's around A. 

Consider now the lower bound on the survival probability. The probabilities Q2 and Q^ are not changed by a 
modification of the initial distribution of traps around the particle. The probability Qi that at t = a box of size 
L around A contains no traps is changed. It is again easiest to start with a discrete lattice and define Con as the 
probability that there is a gap of n sites to the right of A given that there is a single trap to the right of A as 
determined by the distribution qi. Thus 

n n n 

w™ = Yl{l-q;)^\nUn^^H'^-qr)~ -^q^■ (A7) 

i=0 i=0 4=0 

The continuum limit gives 

a}(x) = exp ( - / q{r)dr j . (A8) 

If there are N = pR to the right of A, with i? ^ 00, the probability of a gap of size x to the right of A is 
Lo(x) = uj^ (x) = ujP^{x). Inserting the expression for q{r) and carrying out the integrals for R —> 00 gives 

fexp [-pxV(2(r))] , x < {r) 
uj{x) = I (A9) 

[exp[-p{x- (r)/2)], x> (r) 

The same result is valid for a gap to the left of A, so that the probability Qi for there to be a gap of length L centered 
at A then is 

gi=exp[-piV(4(r))] , for L/2<{r) (AlO) 

Qi = exp [~p{L - (r))] , for L/2 > (r) (All) 

The lower bound on the survival probability then is 

Pm = Qi-exp \^^r- - ^P\l—) ' (A12) 

where we have inserted the explicit expressions for Q2 and Q^ found earlier. As before, one can choose the box size 
L to maximize this lower bound. There are now two possibilities. One is to assume that L < 2(r), whence 

. X 4 / TT^D't AL^ l4:Dt\ , . , 

P4t)^-exp(^-^^-p^-2p^— j. (A13) 

1/4 
Maximizing with respect to L gives L* = (A-n'^D'{r)t/p^ , which grows with time and is therefore in contradiction 

with the assumption L < 2(r). Thus we have to take the other possibility, L> 2(r), whence Eq. (|A11|I leads to 

r. , . 4 / n^D't , , l4:Dt\ , . , 

PL{t) = - exp I j^ -pL + p{r) -2p^—-U (A14) 

and maximizing with respect to L gives L* = (2TT^D't/ p) . For t — > 00, the new term p{r) in the exponent is 
negligible compared to the terms that grow with time. The asymptotic lower bound is therefore identical to that 
obtained by Bray and Blythe for the Poisson distribution of traps around A. 

Although the distribution q{r) is not known in the subdiffusive case, it will surely be the case that at large distances 
from A it is independent of r even if the trapping process has been going on for a long time r before the start of the 
observations at f = 0. Therefore in this case as well the asymptotic survival probability results obtained in this paper 
are applicable. 
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